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ABSTRACT 


A  vlscous-lnvlscld  Interaction  method  for  three-dimensional  flows,  in 
which  the  partla lly-parabolic  Reynolds  equations  are  coupled  with  an  Invlscld- 
flow  solution  procedure  in  an  Interactive  and  iterative  manner,  is  applied  to 
two  simple  three-dimensional  bodies  for  which  experimental  data  are  available 
for  comparison.  The  relative  merits  of  interactive  and  global  solution  pro¬ 
cedures  are  evaluated  by  comparing  the  vlscous-lnvlscld  interaction  solutions 
with  noiilnteractlve  large-domain  solutions  of  only  the  viscous-flow  equa¬ 


tions.  Both  methods  yield  satisfactory  results,  although  the  interaction 
solutions  are  found  to  be  computationally  more  efficient  for  the  cases  con¬ 
sidered. 
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I.  INTRODUCTION 


Recently,  the  authors  (Stern  et  al.,  1986)  described  a  viscous-lnvlscld 
interaction  method  for  three-dimensional  flows  which  require  equations  more 
general  than  the  classical  boundary-layer  equations.  The  method  couples  the 
partially-parabolic  Reynolds  equations  with  an  Invlscld-f low  solution  pro¬ 
cedure  in  an  interactive  and  Iterative  manner.  The  relative  merits  of  inter¬ 
active  and  global  solution  procedures  were  evaluated  by  comparing  the  Interac¬ 
tion  solutions  and  large-domain  solutions  of  only  the  viscous-flow  equations 
with  experimental  data  and  other  computational  methods  for  two-dimensional  and 
axisymmetrlc  bodies.  Both  methods  were  found  to  be  satisfactory,  and  the 
results  clearly  demonstrated  the  feasibility  of  higher-order  viscous-lnvlscld 
Interaction  procedures. 

This  report  describes  the  results  from  applications  of  the  aforementioned 
methods  to  two  simple  three-dimensional  bodies.  Since  the  computational 
method  is  described  In  the  previous  report  It  suffices  here  to  provide  an 
outline.  However,  during  the  course  of  the  present  calculations  It  was  found 
that  modifications  of  some  of  the  finite-difference  formulas  were  required  to 
deal  with  three-dimensional  bodies.  These  are  provided  here.  Also,  correc¬ 
tions  to  errors  in  the  text  of  Stern  et  al.  (1986),  discovered  after  printing, 
are  given  in  Appendix  T. 

II.  OUTLINE  OF  THE  CIWPDTATIONAL  METHOD 

fionsider  the  flow  past  a  streamlined  three-dimensional  body  fixe!  in  a 
uniform  stream  of  velocity  U^.  The  classical  approach  to  the  analysis  of  such 
a  flow  field  is  to  divide  it  into  three  regions,  as  shown  In  figure  1,  In  each 
of  which  a  different  set  of  approximations  Is  made  In  the  Navler-Stokes  equa¬ 
tions  for  laminar  flow,  and  the  corresponding  Reynolds  equations  for  turbulent 
flow.  The  resulting  equations  for  the  tnvlscld  region  1  and  the  boundary- 
laver  region  2  are  well  known.  In  the  absence  of  flow  "eversal  at  the  tall, 
the  thick  stern-flow  boundary  layer  and  the  wake  (region  3)  can  be  described 
bv  rhi>  so-called  par  1 1  a  1  1  y-parabo  1  i  c  equations  in  which  only  streamwlse  gradi¬ 
ents  of  the  viscous  and  turbulent  stresses  are  neglected.  The  interaction  (or 
lack  of  one)  between  regions  1  and  ?  is  the  classical  problem  of  boundary- 
layer  theory.  In  the  present  work,  we  are  concerned  with  the  flow  In  region  3 
and  its  interaction  with  that  in  region  1. 
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There  are  two  possible  approaches  to  the  solution  of  the  partially-para- 
bolic  equations  for  region  3:  a  global  approach  in  which  one  set  of  governing 
equations  appropriate  for  both  the  Inviscid-  and  viscous-flow  regions  is 
solved  using  a  large  solution  domain  (regions  3  and  1)  so  as  to  capture  the 
entire  zone  of  viscous-invlscld  interaction;  and  an  Interactive  approach  in 
which  different  sets  of  governing  equations  are  used  for  each  region  and  the 
complete  solution  obtained  through  the  use  of  an  interaction  law,  i.e.  ,  patch¬ 
ing  or  matching  conditions.  The  former  approach  is  somewhat  more  rigorous 
since  it  does  not  rely  on  the  patching  or  matching  conditions  which  usually 
Involve  further  approximations.  Nonetheless,  for  a  variety  of  reasons,  both 
types  of  approaches  are  of  Interest.  In  the  present  work,  both  are  used  in 
order  to  assess  the  accuracy  of  the  interaction  solution  and  highlight  the 
essential  differences  between  the  two  strategies,  including  their  computa¬ 
tional  efficiency. 

Besides  the  size  of  the  solution  domain,  the  primary  difference  between 
the  two  approaches  lies  in  the  conditions  specified  on  the  outer  boundary 
(see  figure  1).  For  the  noninteractlve  large-domain  solutions,  uniform-flow 
conditions  are  appropriate.  The  interactive  solutions  require  the  speci¬ 
fication  of  the  match  boundary  as  well  as  an  interaction  law,  and  also  a 
method  for  calculating  the  inviscid  flow.  In  the  previous  work,  solutions 
were  obtained  with  the  match  boundary  just  beyond  the  thickness  of  the  bound¬ 
ary  layer  6.  In  the  present  work,  the  Influence  of  the  location  of  the  match 
boundary  was  investigated  by  obtaining  solutions  both  with  the  match  boundary 
just  beyond  6  and  with  the  match  boundary  at  about  26.  The  interaction  law  is 
based  on  the  concept  of  displacement  thickness.  The  inviscld-f low  solutions 
are  obtained  using  a  conf ormlng-pane 1  source-distribution  method. 

Tn  the  part  la  1 ly-parahol 1 c  equations,  the  velocity  field  is  elliptic  in 
transverse  planes  and  parabolic  in  the  streamwise  direction  while  the  pressure 
field  Is  fully  elliptic.  Solutions  to  these  equations  can  be  obtained  itera¬ 
tively  by  solving  the  parabolic  equations  that  result  when  the  pressure  field 
is  specified  and  subsequently  updating  the  pressure  field  using  the  results 
from  the  parabolic  solutions.  In  the  present  approach,  a  body-fitted  coordin¬ 
ate  system,  in  which  the  axial  coordinate  is  roughly  aligned  with  the  stream¬ 
lines,  is  used.  The  partlally-parabollc  assumptions  are  made  in  this  prese- 


lected  axial  coordinate  direction.  Both  numerical  and  analytic  grid  gener¬ 
ation  techniques  are  used.  The  governing  equations  are  derived  In  nonortho- 
gonal  curvilinear  coordinates  with  velocity  components  along  the  coordinate 
directions.  The  Reynolds  stresses  are  modeled  using  the  k-e  turbulence 
model.  The  equations  are  discretized  using  finite  differences  and  a  staggered 
grid.  The  SIMPLER  algorithm  is  used  for  velocity-pressure  coupling, 
equations  are  solved  implicitly  by  SOR  by  lines. 

III.  EXTENSIONS  FOR  THREE-DIMENSIONAL  FLOW 

The  methods  described  In  Stern  et  al.  (1986)  are  for  general  three-dimen¬ 
sional  flow.  However,  during  the  course  of  the  present  calculations  it  was 
found  that  modifications  of  some  of  the  finite-difference  formulas  were  re¬ 
quired.  These  are  provided  In  Appendix  IT.  The  equation  numbers  used  here 
correspond  to  those  used  previously  in  Appendix  TIT  of  Stern  et  al.  (1986). 
Referring  to  these  Appendices,  we  note  that  the  modifications  Involve  moving 
certain  terras  from  the  source  term  to  the  main  diagonal  of  the  tridiagonal 
matrix  formed  when  solving  the  y-  and  z-momentum  equations  and  the  turbulence 
model  equations  by  SOR  by  lines. 

Another  departure  from  the  previous  procedures  relates  to  the  determln- 

* 

atlon  of  the  displacement  thickness  6  (x,z).  An  approximate  method  for 

* 

determining  5  for  three-dimensional  bodies  In  terms  of  a  local  flux  balance 
was  outlined  earlier.  For  the  present  applications  to  the  simple  elliptical 
cross-section  bodies,  an  even  more  approximate  procedure  was  used;  that  Is, 
equation  (V-11)  of  Stern  et  al.  (1986)  was  used  to  determine  the  displacement 
thicknesses  along  only  the  major  and  minor  axes,  and  these  were  connected 
under  the  assumption  that  the  displacement  body  is  also  elliptic  in  cross 
section. 

IV.  APPLICATIONS  TO  THREE-DIMENSIONAL  BODIES 

A.  Elliptical  Cross-Section  Bodies 

The  two  calculation  schemes  were  applied  to  two  simple  three-dimensional 
bodies:  namely,  the  DTNSRDC  2:1  and  3:1  elliptical  cross-section  bodies  for 
which  experiment il  data  have  been  obtained  by  Huang  and  associates  (Groves  et 
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al. ,  1982;  Huang  et  al. ,  1983  and  1985).  Figure  2  shows  a  schematic  view  of 
the  bodies.  These  bodies  have  the  same  length  and  longitudinal  distribution 
of  cross-sectional  area  as  the  Afterbody  1  for  which  results  were  presented  In 
the  previous  report. 

B.  Calculations 

The  calculations  for  both  bodies  were  performed  for  a  body-length  Rey¬ 
nolds  number  Rn  =  6.5  x  10^  which  corresponds  to  the  experimental  condi¬ 
tions.  Typical  large-  and  sraall-domaln  grids  used  in  the  calculations  are 
shown  in  figures  3  and  4  for  the  2:1  and  3:1  bodies,  respectively.  Figures 
3a, b  and  4a, b  show  the  grids  in  the  minor  axis  plane  while  figures  3c, d  and 
4c, d  show  the  cross-sectional  views  at  X  =  .8966.  For  the  2:1  body,  the 
large-domain  grid  was  obtained  using  a  numerical  technique  through  the  solu¬ 
tion  of  elliptic  partial  differential  equations  whereas  the  small-domain  grid 
was  constructed  using  an  analytic  technique  in  which  grids  for  each  cross 
section  were  pieced  together  (stacke'^'  In  the  streamwise  direction.  For  the 
3:1  body,  both  the  large-  and  small-domain  grids  were  obtained  using  the 
numerical  technique.  In  this  case,  the  difference  between  the  two  grids  Is 
chat  Che  small-domain  grid  was  obtained  simply  by  deleting  that  portion  of  the 
large-domain  grid  which  lay  beyond  about  two  boundary-layer  thicknesses.  This 
was  done  to  evaluate  the  accuracy  of  the  two  approaches  and  to  determine  the 
effects  of  the  placement  of  the  outer  boundary  in  the  interaction  solution. 

Referring  to  figures  1  and  3a, c  for  notation,  for  the  2:1  body  large- 
domain  solutions,  there  are  55,  18,  and  5  grid  points  in  the  axial,  radial, 
and  glrthwlse  directions,  respectively.  The  solution  domain  is  bounded  by  the 
inlet  at  =  .5,  the  exit  Sp-  at  x^  =  2.5,  the  outer  boundary  at  y^  =  1.05, 
the  body  surface,  and  the  two  planes  of  symmetry  (0  =  0“  and  90°).  The  first 
grid  point  off  Che  body  surface  was  located  in  the  range  150  <  y"*"  <  300. 
Referring  to  fl)Vires  1  and  3b, d,  for  the  small-domain  solutions,  there  are  55, 
15,  and  5  grid  points.  The  inlet  and  exit  boundaries  as  well  as  the  first 
grid  point  off  the  body  surface  have  the  same  values  as  in  the  large-domain 
grid,  but  the  outer  boundary  is  now  placed  at  y^  =  where  6(x,z)  is  the 
boundary  layer  thickness  which  was  specified  from  the  experimental  results. 
In  both  cases,  the  grid  expansion  was  specified  such  that  the  first  two  grid 
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points  He  in  the  log-law  region.  For  the  3:1  body  (figures  1  and  4),  the 
number  of  grid  points  and  boundary  values.  Including  the  location  of  the  first 
grid  point  off  the  body  surface,  are  similar  to  those  In  the  the  large-domain 
grid  for  the  2:1  body.  The  difference  between  the  large-  (figures  4a, c)  and 
small-domain  (figures  4b, d)  grids  for  the  3:1  body  is  that  In  the  former  19 
radial  grid  points  were  used  whereas  In  the  latter  only  11  were  used.  This 
places  the  outer  boundary  in  the  Interaction  solution  at  about  twice  the 
boundary-layer  thickness. 


For  both  applications  standard  turbulent  flat-plate  profiles  were  used  to 
specify  the  initial  conditions.  All  other  boundary  conditions  were  prescribed 
as  discussed  In  Section  TV.D  of  Stern  et  al.  (1986).  The  calculations  were 


performed  on  a  Prime  9950  minicomputer.  An  underrelaxation  factor  was  used 


In  solving  the  pressure  equation.  For  the  large-domain  solutions  ot  =  .2-. 4 

P 


whereas  for  the  Interaction  solutions  a  =  .3-. 4.  The  large-domain  solution 

P 


for  the  2:1  body  converged  In  40  iterations  and  required  120  minutes  of  cpu 
time,  and  for  the  3:1  body,  it  converged  In  45  Iterations  and  required  135 
minutes  of  cpu  time.  The  Interaction  solutions  converged  in  30  iterations  and 
required  90  minutes  of  cpu  time  In  both  cases. 


C.  Results  and  Discussion 


The  results  for  the  2:1  body  are  presented  in  figures  5  through  13  and 
those  for  the  3:1  body  In  figures  14  through  22.  First,  the  convergence 
history  of  both  approaches  is  considered.  Figures  5  and  6  show  the  conver¬ 
gence  histories  with  regard  to  the  pressure  distribution  for  the  2:1  body  for 
the  large-domain  and  interaction  solutions,  respectively.  For  the  former, 
values  are  shown  for  every  10  global  iterations,  whereas  for  the  latter, 
values  are  shown  after  every  5  global  Iterations.  Figure  7  shows  the  conver¬ 


gence  history  of  the  displacement  thickness  for  the  interaction  solution. 


Corresponding  results  for  the  3:1  body  are  shown  in  figures  14  through  16. 
The  interaction  solutions  were  started  with  free-stream  edge  conditions  (H^  = 
1.,  =  0.,  Pg  =  0).  After  15  1;  atlons  for  the  2:1  body  and  20  iterations 
for  the  3:1  body,  the  edge  conditions  were  updated  using  the  latest  displace¬ 


ment  thickness.  Subsequent  updates  of  the  edge  conditions  were  not  made. 


Note  that  for  Afterbody  1  it  was  shown  that  small  changes  In  displacement  body 


P5; 
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after  intermediate  convergence  was  achieved  (circa.  15-20  iterations),  led  to 
insignificant  changes  in  the  edge  conditions.  This  was  also  confirmed  for  the 
present  applications  through  numerical  experimentation.  The  magnitude  of  the 
edge  velocity  Qg  =  /l-^  used  in  the  interaction  solutions  is  shown  in  figures 
8  and  17  for  the  2:1  ind  3:1  bodies,  respectively.  As  was  the  case  for 
axlsymmetric  flow  and  discussed  with  reference  to  figures  27  and  28  of  Stern 
et  al.  (1986),  the  two  solutions  show  quite  different  convergence  character¬ 
istics.  The  large-domain  solutions  converge  monotonically  in  about  40  to  45 
Iterations  for  both  bodies,  whereas  the  Interaction  solutions  converge  with 
oscillations  in  about  30  iterations  and  with  two  stages  (for  both  cases).  The 
first  stage  is  with  free-stream  edge  conditions  and  leads  to  an  underpredlc- 

if 

tion  of  both  and  S  .  The  second  stage  Is  with  the  displacement-body  edge 
conditions  and  the  solutions  converge  quite  rapidly.  Figures  9  and  18  show 
the  final  values  of  displacement  thickness  for  the  interaction  solutions 
For  =  0°  and  90°  for  the  2:1  and  3:1  bodies,  respectively.  Note  that,  in 
both  cases,  the  displacement  thickness  is  larger  along  6  =  90°  than  along  9  = 
0°  leading  to  stronger  viscous-lnvlscid  interaction  along  9  *  90°. 

Next,  a  comparison  is  made  between  the  two  solutions  and  also  the  exper¬ 
imental  data.  Figure  10  shows  the  pressure  distribution  on  the  surface  of  tiie 
body  and  along  the  wake  centerline  for  the  2:1  body  for  three  glrthwlse  posi¬ 
tions  9  =  (0°,  80°,  90°),  0=0°  corresponding  to  the  waterplane  (minor  axis) 
symmetry  plane  and  9  =  90°  to  the  keel  (major  axis)  symmetry  plane  (see  figure 
2).  The  converged  results  from  the  two  methods,  l.e.,  the  noninteractive 
large-domain  solutions  and  the  sma 1 1 -domain  vlscous-lnvl scid  interaction 
solutions,  are  compared  with  the  experimental  data  and  the  inviscid-f low 
solution  without  Interaction.  It  is  seen  that  both  methods  are  in  good  agree¬ 
ment  with  the  experimental  data.  However,  just  as  was  the  case  for  the 
,1  s i s vmmot r 1 c  bodies,  the  interaction  solution  is  In  slightly  better  agree- 
mo  M  .  This  is  due  to  the  influence  of  the  initial  conditions  and  the  better 
grid  resolution  within  the  boundary-layer  region  for  the  interaction  solu¬ 
tion.  The  degree  of  viscous-lnvlscid  interaction  due  to  the  stern  boundary 
l  iver  and  wake  for  this  body  Is  similar  in  magnitude  to  that  for  the  parent 
axisvmmetric  shape  (Afterbody  1).  Also,  note  the  larger  viscous  effects 
along  '  -  art'’  than  along  9  =  0°.  Figure  11  shows  the  glrthwlse  pressure 


distribution  at  several  axial  stations.  Consistent  with  figure  10,  the  inter¬ 
action  solution  shows  better  agreement  with  the  experimental  data.  The 
increase  in  pressure  towards  9  =  90“  is  expected  due  to  the  rapid  change  of 
body  curvature  in  this  region.  The  results  for  the  3:1  body  corresponding  to 
figures  10  and  11  are  shown  in  figures  19  and  20.  It  is  seen  that  the  vls- 
cous-lnvlscld  interaction  is  larger  for  the  3:1  body  than  it  is  for  the  2:1 
body,  especially  near  the  keel.  The  differences  between  the  large-domain  and 
interaction  solutions  for  the  3:1  body  are  less  pronounced  than  those  for  the 
2:1  body  and  previously  for  axisymmetrlc  bodies.  This  is  no  doubt  due  to  the 
fact  that  for  the  3:1  body  both  methods  have  the  same  number  and  distributions 
of  radial  grid  points  within  the  viscous-flow  region.  The  comparison  with  the 
experimental  data  is  not  as  good  as  shown  for  the  2:1  body. 

The  calculated  distributions  of  the  wall-shear  velocity  along  three 
glrthwlse  positions  (9  =  0“,  77“ ,  90“ )  on  the  2;1  body  are  compared  with  the 
experimental  values  deduced  from  Preston  tubes  in  figure  12.  In  this  case, 
the  agreement  between  the  two  sets  of  calculations  is  good  but  the  experimen¬ 
tal  values  are  lower  especially  in  the  region  X  >  0.9.  This  discrepancy  could 
be  attributed  to  deficiencies  of  the  turbulence  model;  however,  as  noted  by 
Chen  and  Patel  (1985)  from  a  reanalysls  of  the  data,  the  values  deduced 
from  the  measured  velocity  profiles  and  Clauser  plots  are  higher  and  are  in 
better  agreement  with  the  calculations.  The  corresponding  results  for  the  3:1 
body  are  shown  in  figure  21.  Consistent  with  the  earlier  discussion,  the  two 
solutions  are  in  even  closer  agreement  than  that  in  the  case  of  the  2:1 
body.  Unfortunately,  experimental  values  are  not  quoted  to  verify  the 
resul ts. 

Figure  13  shows  the  profiles  of  (U,V,W,C  ,k,e)  for  the  2:1  body  at  five 

P 

stations  on  the  body  X  =  (.767,  .894,  .914,  .944,  .956)  where  experimental 
data  is  available  (see  figure  2).  In  these  figures,  the  radial  coordinate  Is 
defined  as  Y  =  (R  -  R„)/R„„,.  where  R_(x)  is  the  local  body  radius,  R  is  the 
distance  from  the  body  surface  along  the  local  normal  to  the  body  cross  sec¬ 
tion,  and  Is  the  maximum  value  of  the  geometrical  mean  radius  (i.e.,  the 
maximum  radius  of  the  Afterbody  1).  The  calculated  velocity  components  have 
been  Interpolated  and  transformed  into  the  experimental  coordinate  system  in 
order  to  make  the  comparisons.  The  U  component  Is  In  the  axial  cartesian 


coordinate  direction,  the  V  component  is  in  the  direction  of  the  local  normal 
to  the  body  cross  section,  and  the  W  component  is  orthogonal  to  both  U  and 
V.  Consistent  with  the  results  discussed  above,  both  solutions  are  in  good 
agreement;  however,  some  differences  are  evident,  especially  in  the  pres¬ 
sure.  That  is,  while  both  solutions  Indicate  similar  transverse  pressure 
gradients  Py  there  are  systematic  differences  in  pressure  magnitude.  The 
large-domain  solution  predicts  lower  pressure  than  the  interaction  solution 
which,  in  general,  shows  better  agreement  with  the  experimental  data.  A 
similar  trend  was  found  previously  in  the  results  for  axlsymmetrlc  bodies. 
This  may  be  related  to  the  outer  boundary  conditions  in  the  large-domain  solu¬ 
tion  since  its  characteristics  are  similar  to  a  blockage  effect.  The  velocity 
components  are  in  good  agreement  with  the  experimental  data.  This  is  espe¬ 
cially  true  for  the  U  component.  The  V  component  is  slightly  over  predicted. 
The  W  component  is  quite  small  and  difficult  to  compare  with  the  experimental 
values.  Actually,  the  experimental  W  values  indicate  an  error  in  the  align¬ 
ment  of  the  model  since  nonzero  values  are  shown  for  the  symmetry  planes,  9  = 
0°  and  90°.  The  calculations  show  a  systematic  overpredict  ion  of  the  turbu¬ 
lent  kinetic  energy  as  the  boundary  layer  thickens;  namely,  as  the  tall  is 
approached  and  Increasingly  so  towards  9  =  90®.  This  is  undoubtedly  due  to 
deficiencies  of  the  turbulence  model. 

Corresponding  results  for  the  3:1  body  are  shown  in  figure  22.  In  this 
case,  the  body  station  values  are  X  =  (.719,  .810,  .894,  .934,  .954)  (see 
figure  2).  Here  again,  the  two  solutions  are  in  even  closer  agreement  than 
that  shown  for  the  2:1  body;  however,  a  difference  in  pressure  magnitude  is 
still  evident.  The  level  of  agreement  with  the  experimental  data  for  the 
symmetry  planes  (i.e.,  9=0°  and  90°)  is  very  good  and  similar  to  that  shown 
previously  for  the  2:1  body.  However,  the  results  for  9  =  80°  do  not  show 
overall  good  agreement  with  the  experimental  data.  The  U  component  shows  a 
discrepancy  that  begins  in  the  outer  part  of  the  boundary  layer  at  X  =  .81  and 
becomes  progressively  worse  and  encompasses  more  of  the  boundary  layer  as  the 
tail  is  approached.  The  V  component  is  overpredicted,  the  W  component  indi¬ 
cates  negative  values  whereas  the  experimental  values  are  positive,  and  the 
turbulent  kinetic  energy  is  overpredicted.  A  similar  discrepancy  in  the  U 
component  occurs  in  the  outer  part  of  the  boundary  layer  in  the  previous 
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results  for  the  2:1  body  and  Afterbody  5;  however,  in  these  cases  It  Is  con¬ 
fined  to  the  outer  part  of  the  boundary  layer  and  the  error  is  small.  Further 
numerical  testing  is  required,  Including  grid  refinement  and  comparisons  with 
more  extensive  sets  of  experimental  data,  in  order  to  determine  the  cause. 

V.  CONCLUDING  REMARKS 

A  viscous-lnvlscld  interaction  method  for  three-dimensional  flows  which 
require  higher-order  viscous-flow  equations,  has  been  validated  by  comparisons 
with  experimental  data  for  two  simple  three-dimensional  bodies;  namely,  the 
DTNSRDC  2:1  and  3:1  elliptical  cross-section  bodies.  The  relative  merits  of 
interactive  and  global  solution  procedures  were  evaluated  by  comparing  the 
interaction  solutions  with  the  large-domain  solutions  of  only  the  viscous-flow 
equations.  Both  methods  yield  similar  results  in  good  general  agreement  with 
the  data.  The  level  of  agreement  between  the  two  approaches  and  with  the 
experimental  data  is  similar  to  that  shown  previously  for  axisyrametric  bodies. 

As  was  anticipated  on  the  basis  of  previous  calculations  for  axisymmetrlc 
bodies,  the  interaction  solution  is  computationally  more  efficient  than  the 
large-domain  solution,  l.e.,  25%  savings  in  overall  cpu  time  for  the  cases 
considered.  This  is  due  to  the  fact  Chat  the  prescription  of  the  Invlscld- 
flow  at  the  boundary-layer  edge  accelerates  the  convergence  of  the  viscous- 
flow  solution.  Although  there  is  an  additional  computational  effort  in  cal¬ 
culating  the  Invlscid  flow,  the  savings  in  the  calculation  of  the  viscous-flow 
is  more  substantial. 

The  feasibility  of  both  the  present  approaches  has  now  been  demonstrated 
for  various  levels  of  flow  field  complexity:  two-dimensional,  axisymmetrlc, 
and  simple  three-dimensional  bodies.  The  results  are  quite  satisfactory  and 
warrant  further  work  in  extending  the  present  procedures  for  more  realistic 
ship  geometries  and  complex  flow  fields.  Of  particular  interest,  is  the 
extension  of  the  Interaction  approach  for  the  calculation  of  ship  boundary 
layers  and  wakes,  including  the  effects  of  free-surfce  waves,  i.e.,  for  non¬ 
zero  Froude  numbers,  by  coupling  the  present  viscous-flow  method  with  a  poten¬ 
tial-flow  soltitlon  procedure,  such  as  Dawson's  method  (1977),  in  a  similar 
manner  as  has  been  done  presently  for  zero  Froude  number.  It  is  believed  that 
an  interaction  approach  Is  more  readily  extendible  to  the  problem  of  calculat¬ 
ing  ship  boundary  layers  at  nonzero  Froude  numbers. 
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I’lqure  1.  Dutinitioii  Sketch  of  Flow-Field  Regions. 
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APPENDIX  I:  Corrections  to  Stern  et  al.  (1986) 

This  appendix  provides  corrections  to  errors  in  the  text  of  Stern  et  al. 
(1986)  which  were  discovered  after  printing. 

1)  p.  10  line  1,  equation  (III-3)  should  read: 

*  ^ 

5  =  f  (1  -  U/U  )  dy 

e 
o 

2)  p.  32  line  5,  equation  (V-6)  should  read: 


G  =  -  1/R 


R  =  X  -  C 


3)  p.  32  line  10,  equation  (V-7)  should  read: 


c  dG 

2na  -  J  o  dS  =  -  U  n, 

^  3n  o  1 

B 

A)  p.  35  line  16,  in  equation  (V-12)  o  is  defined  by: 


1  d  * 

o  -  4-  ^  (U  6  ) 
2Tf  dx  e 


5)  p.  112  line  19,  the  equation  should  read: 


,  ^  3U  9V  3W  , 

(U^,  V^,  W^)  -  g^^) 

2  2 

6)  p.  115  line  7,  in  equation  (B-1)  the  term  R^U  should  read; 


7)  P.  117  lines  18  and  19,  equation  (B-9)  should  read: 

1  *^1 

^12  "  7  ^~2  ^AUbjj+  H(Ubj2+  V2)  +  G(Ubj3+  Vg)} 
s 


«*.  •  .  . 

"  -  .•  .  •.* 


/‘v‘v '.-v^v/S  > 


A’.  .*/. 


•V  •  •' 


V  •" 


■v 


■r.  .v.J-' 


33 


-  ■  ^  K  ' 


+  -j  {H(Uj+  Uajj)  +  B(U2+  Uaj2)  +  F(U3+  Uaj^)} ] 
s 

8)  p.  117  lines  20  and  21,  equation  (B-10)  should  read; 

1  ^1 

^13  "  T  f“2  ^AUCjj+  8(00^2"^  ^2^  G(Ucj3+  W3)} 
s 

+  -^  {GvUj+  Uajj)  +  F(U2+  Uaj2)  +  0(03+  Uaj3)}] 


(B-9) 


y.v^- 


"v-' 


v' 

.V 

5-:;5 

m 


(B-10) 


9)  p.  118  lines  2  and  3,  equation  (B-12)  should  "ead: 


APPENDIX  II.  Modifications  of  Equations  In  Discretized 
Fora  for  Three— Dlaenslonal  Flow 

As  discussed  In  Section  III,  some  modifications  of  the  finite-difference 
formulas  presented  in  Appendix  III  of  Stern  et  al.  (1986)  were  found  to  be 
necessary  upon  performing  calculations  for  three-dimensional  bodies.  These 
are  provided  below; 


1)  p.  126  lines  1  and  2,  in  the  right  hand  side  of  equation  (C-15)  e 
should  be  replaced  by  ^23’  where 

T  =  e  -  -L. 

23  23  2h,  3z 


23 


2)  p.  126  after  line  5,  the  following  equation  should  be  Inserted: 

3 


3  ,  3V  , 

(V.  Z~TZ) 


V 


_  _  ._jt .  3v  ^t  3^v 

h  3z  ^  "t  h  Sz''  h  3z  ^h  3z  2.2 
3  3  3  3  h~  9z 


V  il  -  — 

"  2  ^h '3z  ^h"^^  1  ^\,n+l  ”  \,n-l  ^ 

3  3m-  y.n 


1 


(—1  f  +1"  2  +  V*  ,)  =  yd.  -  yd,/  (C-16-1) 

2  1  m,n+l  m,n  m,n-l  5  ^6  m,n 


m  -  y  ,  n 


3) 


4) 


P- 


P- 


128  line  7,  the  right  hand  side  of  equation  (C-21)  should  read: 

128  line  12,  the  right  hand  side  of  equation  (C-22)  should  read: 


(C-21) 


(C-22) 


3)  p.  129  lines  2-13,  should  be  replaced  by: 


Collecting  all  terras  (C-12)  -  (C-23): 


i  £-1  £  £  £ 

ya,(V  -  V  )  +  yb,  V  ,  -  yb«V  +  yb,V  ,  +  yc,  +  yd, 

1  m,n  m,n  •'  1  in+l,n  ^  2  ra,n  3  m-l,n  1  1 

-^•^2  ^  y^^ra.n  '  y^4Cl.n“  ^  ^  ^  yf2'^.n 

-  y^3  Cl.n  ■•  -  y^5  -  ‘  =  0 


(yb3  -  ^*^4  "  y^3^  Vl.n  ^y^2  "  ^^^2  ^^3  ^*^6  •"  ^^2  ^^2^  \.n 

+  (ybj-  yd2-  yfj)  ya^V^"^-  yc^-yd^^  yd,.-  ye^+  yf^+  yf5+  yf^-  ygj+  yh^ 

Finally,  the  coefficients  in  equation  (TV-19)  are: 


bj  =  yb^  -  yd^  -  yf^ 


=  ya^  -  yb^  +  ydj  +  yd^  +  yf^  +  yf^  +  yg. 


(C-24) 


b^  =  ybj  -  yd^  -  yf^ 


Sv  =  yaj  yCj-  ydj+  yd^+  ydj+  yf^+  yf^+  yf^-  ygj+  yh^ 


Pv  =  -  ye 


6)  p.  131  lines  6  and  7,  in  the  right  hand  side  of  equation  (C-28)  e 
should  be  replaced  by  ^23’  where 

7  -  e  -_!-9" 

23  23  2h.  3y 


KSSi 


to 


exi 


“•y^' 


7)  p.  132  after  line  1,  the  following  equation  should  be  Inserted: 

-VI  V  V  ? 

9  ,_t  ^  ^  _3 _  _t  3W  _t  3^ 

^  ^2  h^  3y2 


C'.v.'! 


=  zaj-  zc^-+  zd^+  zd^+  zf^+  zf^+  zg^ 


C3  ■  "r  "'‘2"  "'4 


(C-37) 


za ,  W 
1  ni,n 


zbj-  zdj+  zd^+  zfj+  zf^+  zh^ 


11)  p.  137  line  4,  the  right  hand  side  of  equation  (C-44)  should  read: 


=  Fd,  -  Fd,  F  (CD) 
6  7  m,n 


12)  p.  137  lines  13-15,  should  be  replaced  by: 


I  I  in 

+  Fd_F  -  Fd-F  ,  -  Fd,-  Fd--  Fd,+  Fd,F  -  Fe,  +  Fe.F  =  0 

2m, n  3m-l,n  4  5  6  7m, n  1  2m, n 


-(Fb,+  Fd,)  F*^  ,  +  (Fa,+  Fd,  +  Fd-,+  Fe_)  F^  +  (Fb,-  Fd.)  F^, 

1  3  m-l,n  1  2  7  2  m,n  1  1  nH-l,n 


(C-44) 


13)  p.  138  line  3,  should  be  replaced  by: 


(d,e)2  =  Faj+  Fd2+  Fd^+  Fe^ 


(C-47) 
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